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Abstract 

An explicit Bargmann symmetry constraint is computed and its associated binary non- 
linearization of Lax pairs is carried out for the super Dirac systems. Under the obtained 
symmetry constraint, the n-th flow of the super Dirac hierarchy is decomposed into two 
super finite-dimensional integrable Hamiltonian systems, defined over the supersymmetry 
manifold /j47V|2W -^[^Yi the corresponding dynamical variables x and t„. The integrals of 
motion required for Liouville integrability are explicitly given. 
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1 Introduction 

For almost twenty years, much attention has been paid to the construction of finite-dimensional 
integrable systems from soliton equations by using symmetry constraints. Either (2-|-l)-dimensional 
soliton equations [U El |3] or (l-l-l)-dimensional soliton equations ^ \5\ can be decomposed into 
'Corresponding author, E-mail address:jshe@ustc. edu.cn, hejingsong@nbu.edu.cn 
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compatible finite-dimensional integrable systems. It is known that a crucial idea in carrying out 
symmetry constraints is the nonlinearization of Lax pairs for soliton hierarchies, and symmetry 
constraints give relations of potentials with eigenfunctions and adjoint eigenfunctions of Lax 
pairs so that solutions to soliton equations can be obtained by solving Jacobi inversion problems 
[6] . The nonlinearization of Lax pairs is classified into mono- nonlinearization [3 [HI [9] and binary 
nonlinearization jlO^ \TT\ [T2] . 

The technique of nonlinearization has been successfully applied to many well-known (1-1-1)- 
dimensional soliton equations, such as the AKNS system [1][8], the KdV equation [5] and the 
Dirac system |13) . But there are few results on nonlinearization of super integrable systems, 
existing in the literature. Studies provides many examples of supersymmetry integrable systems, 
with super dependent variables and/or super independent variables [Ml [HI [Ml [171 [TBI [19] . Very 
recently, nonlinearization was made for the super AKNS system [20j and the corresponding super 
finite-dimensional Hamiltonian systems were generated. In this paper, we would like to analyze 
binary nonlinearization for the super Dirac systems under a Bargmann symmetry constraint. 

The paper is organized as follows. In the next section, we will recall the super Dirac soliton 
hierarchy and its super Hamiltonian structure. Then in section 3, we compute a Bargmann 
symmetry constraint for the potential of the super Dirac hierarchy. In section 4, we apply 
binary nonlinearization to the super Dirac hierarchy, and then obtain super finite-dimensional 
integrable Hamiltonian systems on the supersymmetry manifold i?^^'^^, whose integrals of 
motion are explicitly given. Some conclusions and remarks are listed in section 5. 



2 The Super Dirac Hierarchy 



h = U(t>, u 



X + s a \ 



-X + s -r P 
\ P -a / 



(1) 



The super Dirac spectral problem associated with the Lie super-algebra -6(0, 1) is given by [21j 

/ r \ 

s 
a 

V / 

where A is a spectral parameter, r and s are even variables, and a and [3 are odd variables. 
Taking 



V '^3 / 



V 



( C A + B p\ 

A-B -C 6 
\ 6 -p J 



2 



the co-adjoint equation associated with ([I]) 14 = [f/, V] gives 



If we set 



= -2XC + 2rB -ap + (35, 
= 2rA - 2sC - ap - p5, 
= 2\A - 2sB + a<5 + (3p, 
p^ = -(3{A + B)-aC + {\ + s)5 + rp, 
^ 4 = (-A + s)p -r6- a{A - B) + f3C. 



i>0 i>0 

equation ([2]) is equivalent to 



i>0 



i>0 



Ao = Co = po = 5o = 0, 

Ai+i = \Ci^x + sBi - ^a5i - \[3pi, i > 0, 

Ci+i = -\Ai^x + rBi-\api + \f)5i, i > 0, 

Pi+i = -5i^x - r6i + spi - a{Ai - Bi) + /3Ci, i > 0, 

6i+i = pi^x - rpi - s6i + P{Ai + Bi) + aCi, i > 0, 

Bi+i^x = 2rAi+i - 2sCi+i - api+i - (36i+i, i > 0, 



which results in recurrence relations 



where 



C 



(Ci+i,ylj+i,5i+i, -pi+i)'^ = C{Ci,Ai,6i, -piY , i > 0, 
Bi = d-H2rA,-2sCi-api- p5i), i>0, 

( -2rd-^s -\d + 2rd-^r \l3-rd~^l5 ^a + rd'^a 

\d-2sd-^s 2sd-^r -^a-sd'^p + sd-^a 

a-2(3d-h (3 + 2f3d-^r -s-pd'^p -d + r + pd'^a 

\ —P + 2ad^^s a — 2ad^^r d + r + ad~^P s — ad~^a 



Upon choosing the initial conditions 



Ao = Co = po = So = 0, Bo = 1, 



all other Ai, Bi,Ci, pi,5i, i > 1, can be worked out uniquely by the recurrence relations ([5]). The 
first few results are as follows 

Ai = s, Bi = 0, Ci = r, pi = a, 6i = f3, 

A2 = ^r^, B2 = ^{r^ + s^) + a(3, C2 = -^^x, P2 = -Px, h = a^, 
A3 = -^Srcx + ^{r'^ + s'^)s + saP - ^aar, + 

-B3 = -^{rsx - r^s) + aax + PPx, 
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C3 = --r^x + -{r'^ + s'^)r + ra(3 + -aPx - -axP, 

P3 = -axx + ^{r'^ + s'^)a - ^r^a -^SxP - rax - sPx, 
h = -Pxx + ^{r'^ + s'^)P + ^rxP -^Sxa- sax + r/3x. 
Let us associate the spectral problem ([T|) with the following auxiliary spectral problem 

= V^^"V = (A"F)+</), (6) 

with 



/ d A, + B, Pi \ 



-c. 



\ 



Ai - Bi 

-Pi / 

where the plus symbol "+" denotes taking the non-negative part in the power of A. 

The compatible conditions of the spectral problem ([T]) and the auxiliary spectral problem 
([6]) are 

Ut^-Vt^ + [UM''\=^, n>0, 
which infer the super Dirac soliton hierarchy 

Ut„ = Kn = {2An+l,-2Cn+l,6n+l,-pn+lf , n > 0. 

Here nt„ = Kn in ([8]) is called the n-th Dirac flow of the hierarchy. 
Using the super trace identity [HI [22] 



(7) 
(8) 



(9) 



where Str means the super trace, we can have 



Ai+i 
\ -Pi+i ) 



— i/j, Hi 
du 



Bi+2 



dx, i > 0. 



(10) 



Therefore, the super Dirac soHton hierarchy ([8]) can be written as the following super Hamilto- 
nian form: 



where 



J 



( 

-2 


V 



J 



5H„ 



5u ' 

2 \ 
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fill 



is a supersymplectic operator, and is given by (fTUj) . 

The first non-trivial nonlinear equation of hierarchy (|11|) is given by the second Dirac flow 



^2 ■ 



^ -Pxx + ^{r"^ + s^)/? + rPx - sax + \rxf3 - \sxa, 
Oxx - ^(»^^ + + rax + sf3x + \rxa + is^;/?, 

which possesses a Lax pair of U defined in ([1]) and V^'^'^ defined by 

r\ - \sx \^ + s\ + \rx + i(r2 + s^) + a/3 aA - \ 



(12) 



V 



-A2 + sA + ir^-i(r2 + s2)-a/3 

/3A + 



-aA + f3x 



f3\ + 




Remark 1 We consider all differential equations in the real field and explore the Liouville 
integrability on real symplectic manifolds. We did not see any equivalence between the real 
Dirac soliton hierarchy and the real AKNS soliton hierarchy. For example, it is clear that the 
AKNS system of nonlinear Schrodinger equations and the Dirac system of nonlinear Schrodinger 
equations can not be transformed into each other by any real linear transformations. There 
is a similar situation between the super Dirac soliton hierarchy and the super AKNS soliton 
hierarchy, and between the Liouville integrable constrained flows associated with the two super 
soliton hierarchies. 
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3 The Bargmann symmetry constraint 



In order to compute a Bargmann symmetry constraint, we consider the following adjoint spectral 
problem of the spectral problem ([1]): 

/ -r 



X-s 13 \ 



\ 



-X-s 
—a 



/ ^1 \ 



1p2 
V ^3 / 



(13) 



r —a 
-P / 

where St means the super transposition. The following result is a general formula for the 
variational derivative with respect to the potential u (see for the classical case). 

Lemma 1 Let U{u,X) be an even matrix of order m + n depending on u,Ux,Uxxr ' ' o.nd a 
parameter X. Suppose that (j) = {4>e,(j)o)'^ and ip = (V'ejV'o)^ satisfy the spectral problem and the 
adjoint spectral problem 

(t)x = U (u, X)(p, = -U^\u, X)ip, 

where = ((/)i , • • • , (pm) andipe = {ipi, • • • ,ip m) are even eigenf unctions, and (pg — (^(pirn+i^ ■ ' ' i (pm+n) 
andipo = {ipm+ii ■ ■ ■ ,ipm+n) o,re odd eigenf unctions. Then, the variational derivative of the spec- 
tral parameter X with respect to the potential u is given by 

5X (^e,(-l)^(-Vo)(^) 



where we denote 



p{v) 



0, V is an even variable, 

1, V is an odd variable. 



(14) 



(15) 



By Lemma [U it is not difficult to find that 

/ IplCpl - 1p2(p2 \ 



5X 
5u 



(16) 



1pl(p2 + -02 (Al 
1pl(p3 + i^3fp2 
V ^2(p3 - i'3(pl j 

When zero boundary conditions lim|^|^oo ^ = lim|2.|^oo tp = are imposed, we can obtain a 
characteristic property - a recurrence relation for the variational derivative of A: 



C— = X— 

6u 6u ' 



(17) 



where C and ^ are given by (0) and respectively. 



6 



Let us now discuss two spatial and temporal systems: 





( ] 




( \ 


^ r \j + s a \ 


( \ 










= -Xj + s -r (3 




< 




X 


\ 4>3j ) 


\ ^ -a J 





( <Pij \ 



(18) 



( 



( ^Ij \ 



-Xj - s 



-a 



r —a 
-13 J 



Tp2j 
\ ^3i / 



n 

Y^{A,-B,)XY 

n 

\ i=0 



n n 

- E c^x^;-' E 

1=0 j=0 
n 

- E P^>^^ 

i=0 



n—i 



-{V^^))St(u,Xj) 



E C^X]- 
i=0 



( V'l, \ 



.^{A,-B,)X]-' ES^X]- 

i=0 i=0 



n 

- E /^A" 



1=0 



E c^A"-* 

j=0 
n 

- E ^^>^V 



■ E P.A^ 

j=0 





n—i 



I 



( ^Ij \ 



(19) 



where 1 < j < N and Ai, • • • , Atv are N distinct spectral parameters. Now for the systems (jlSp 
and ()19p . we have the following symmetry constraints 

/I A • 

(20) 



6 v-^ ^A, 



6u ^-^ 6u 

i=i 

The symmetry constraint in the case of /c = is called a Bargmann symmetry constraint |12j . 
It leads to an explicit expression for the potential u, i.e., 

r =< ^'i,$i > - < ^-2,^2 >, 
s =< ^'l,^>2 > + < ^'2,^1 >, 
a = - < ^'2, ^3 > + < ^3, ^1 >, 
/? =< ^-1, $3 > + < ^'3.^2 >, 



(21) 



where we use the following notation, 



= {<Pii,--- ,4'iN)'^, = (V'ii,--- ,'4'iN)'^, i = 1,2,3, 
and < •, • > denotes the standard inner product of the Euclidian space R'^ . 



4 Binary nonlinearization 

In this section, we want to perform binary nonlinearization for the Lax pairs and adjoint Lax 
pairs of the super Dirac hierarchy (jlip . To this end, let us substitute (|2ip into the Lax pairs 
and adjoint Lax pairs (jlSp and ()19p . and then we obtain the following nonlinearized Lax pairs 
and adjoint Lax pairs 













/ 




+ s d ^ 




( ] 
















X, + s - 


-f P 












\ <P^j ) 


X 


\ 4>3j ) 




V 


P 


-d / 








< 


( ^Ij \ 






( ^IJ \ 




( ~" 


Xj - s 


P ) 




( V^i. \ 




tp2o 








-Xj - s 


r 




d 










\ j 


X 








\ —a 




J 




V ^3, J 



(22) 



( i'lj \ 



y(")(u,Aj) 
/ 



<p2j 



n 

E C^X^~' 
n 

i=0 



i=0 



i=0 



E aA"-' E ^iA"- 



1=0 



j=0 



E'^.A 

i=0 



n—i 
3 



E PiA 

j=0 



n—i 

'j 



( ^IJ \ 
\ 4>3j ) 



-{V'^^))St{u,Xj) 



E C^X^~ 
j=0 



E(^i + ^^)Af 

i=Q 

n 

- E P^>^^ 
i=0 



Y.{A,-B,)XJ 

n 

E c-^A"-^ 

i=0 
n 

- E ^.A"-* 

i=0 



. \ 

i=0 
n 




( ^ij \ 


- E Pi>^r 

i=0 




Tp2j 







[ ^3. / 



(23) 
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where 1 < j < and P means an expression of P{u) under the exphcit constraint (j2ip . Note 
that the spatial part of the nonhnearized system ()22p is a system of ordinary differential equations 
with an independent variable x, but for a given n (n > 2), the t„-part of the nonhnearized system 
()23p is a system of ordinary differential equations. Obviously, the system ()22p can be written as 

^l,x = (< ^1, ^-1 > - < ^'2, ^2 >)^1 + (A+ < ^'l, ^>2 > + < 'I'2, ^1 >)^'2 + (- < ^2, ^3 > 
+ < 1-3, $1 >)$3, 

^2,x = (-A+ < ^1, $2 > + < ^2, ^'l >)^'l - (< ^'l, $1 > - < ^2, ^2 >)$2 + (< ^1, ^"3 > 
+ < ^'3,$2 >)^3, 

^3,x = (< 1-1, «>3 > + < ^3, ^2 >)^1 - (- < 1'2, $3 > + < 1'3, ^1 >)$2, 

< 

= -(< ^1,^1 > - < ^2,^'2 >)^'l + (A- < ^'i,$2 > " < ^2,^1 >)^1 + (< 'I'l,$3 > 
+ < 1'3,$2 >)^'3, 

*2,x = -(A+ < ^1, ^-2 > + < 'I'2, ^1 >)^1 + (< 'I'l, $1 > - < 'I'2, ^2 >)*2 - (- < ^'2, $3 > 
+ < 1'3,$1 >)^'3, 

*3,x = -(- < ^2,^3 > + < ^3,«'l >)*1 - (< ^1,«'3 > + < ^3,«'2 >)^'2, 

(24) 

where A = diag{Xi, ■ ■ ■ ,Xn)- When n = 1, the system ([231) is exactly the system (f22]) with 
ti = x. When n = 2, the system (f23]l is 

' ^>i,t2 = (rA - + (A^ + sA + if^ + i(f2 + s^) + d/3)$2 + (aA - /3^.)^3, 

$2,i2 = (-A2 + sA + if^ - + - - (f A - is^)$2 + (/5A + d^)$3, 

$3,t2 = (/3A + d:,)$i - (dA-/3^.)^2, ,,,, 
< _ _ (25) 

^i^t, = -(f A - is^)^'! + (A2 - sA - if^ + i(f2 + s2) + + {ISA + d^)^'3, 

^-a.ta = -(A2 + sA + if^. + i(f=2 + §2) + d/3)^'i + (f A - is^.)^2 - (aA - 4)^'3, 

, ^3,t2 = -(aA - /3^)^'i - (/3A + d^)^'2, 

where f,s,a,P denotes the functions r,s,a,P defined by the explicit constraint (pT]) . and 
fx,Sx, ax,Px are given by 

= 2 < A"!-!, ^>2 > +2 < A^'2, $1 > +2 < 'fi, $2 >^ -2 < ^'2, ^1 
s^ = -2< A^-i, > +2 < A^2, ^"2 > +2(< ^1, $1 > - < ^'2, <^2 >)(< ^2, > - < 'I'l, $2 >), 

< 

=< A^-i, $3 > + < A^-s, ^2 > +(< 'I'l, $2 > - < ^2, ^"1 >)(< ^1, ^>3 > + < ^3, ^2 >), 

^ Pec =< A^2, ^'3 > - < A^3, ^1 > -(< ^1, $2 > - < ^2, ^1 >)(- < ^-2, ^"3 > + < *3, ^^l >), 

which are computed through using the spatial constrained flow (j24p . 

In what follows, we want to prove that the system (j22p is a completely integrable Hamiltonian 
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system in the Liouville sense. Furthermore, we shall prove that the system (j23p is also completely 
integrable under the control of the system (j22p . 

On the one hand, the system (j22p or (j24p can be represented as the following super Hamil- 
tonian form 



dHi ^ dHi ^ dHi 



dHi 



dHi ^ dHi 



5^2 



9$. 



where 



(26) 



1 



1 



< A^i, ^>2 > - < A*2, ^1 > +2(< ^1 > - < ^2, ^2 >) +-(< ^-i, $2 > + < *2, ^1 

+ (- < ^2,^3 > + < ^3,^1 >)(< ^1,^'3 > + < ^3,^2 >)■ 



In addition, the characteristic property (jl7p and the recurrence relations ([5]) ensure that 

ij+i =< A*^'i,^'2 > + < A*^'2,^i >, « > 0, 
Bi+i =< A'^2, > - < A'^-i, $2 >, i > 0, 
Ci+i =< A*^'i,^>i > - < A*^'2,$2 >, ^ > 0, 
5i+i =< A'^i, ^>3 > + < A'^3, ^>2 >, i > 0, 
^ Pi+i = - < A*^'2, ^3> + < A*^3, ^'i >, i> 0. 



(27) 



Then the co-adjoint representation equation Vx = [U, V] remains true. Furthermore, we know 
that the equality = [U, V^] is also true. Let 
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(28) 



Then it is easy to find that = 0, that is to say, F is a generating function of integrals of 
motion for the system ([22|l or ([24ll . Due to F = X^„>o-fnA^", we obtain the following formulas 
of integrals of motion 

-^0 = ~2^0^ -^1 = —BqBi, Fn = —BoBn + - ''^{AiAn-i — BiBn-i+CiCn-i+'^Pi^n-i), n > 2. 

i=l 

(29) 

Substituting ()27p into the above formulas of integrals of motion, we obtain the following expres- 
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sions of Fm{m > 0): 

Fo = Fi =< ^-1,^2 > - < *2,«'l >, 

n-l 

1=1 

+i(< A^-i^i, $1 > - < ^2)(< A""^"^^i, $1 > - < A"-'-i^2, ^2 >) 

+ (- < A'-^<f2,^3 > + < A*~^^'3,«'l >)(< A"-*-l^i,$3 > + < A"-*-1^3,$2)], 

n > 2. (30) 

On the other hand, let us consider the temporal part of the nonlinearized system (j23p . Mak- 
ing use of ([27|1 and (|30]) . the system (f23ll can be represented as the following super Hamiltonian 
form 

_ dFn+i _ dFn+i _ dFn+i 

^1. ^2. ^3. =^ 

This can be checked pretty easily. For example, we can show the last but one equality in the 
above system as follows: 



^2,t„ = -X;(ii + S,)A"-^M/i+^C'iA"-^M'2-5^M""'^3 

1=0 i=0 i=0 

n n 

= -A"^i - 2 J] < A^-^^'2,$i > A"-*^'! + ^(< A*-i^i,$i > - < A'-^^2,^2 >)A" 



i=l 1=1 

3 



+ J](< A'-^^2,^3 > - < A*-i^'3,$i >)A"-^^' 

i=l 

In order to further show the Liouville integrability for the constrained flows (j22p and (j23p . 
we need to prove the commutative property of the integrals of motion {-Fm,}m>Oi under the 
corresponding Poisson bracket: 

X — i J — i 

At this time, we still have an equality t4„ = [^^""^^J) and after a similar discussion, we know 
that F is also a generating function of integrals of motion for (123p . Hence F^i, m > 0, are 
integrals of motion for the system (j23p or (I3ip , which implies that 
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d 

{Fm+i,F„+i} = — F^+i = 0, m,n>0. (33) 
The above equality (p3]) shows that {Fm}m>o are in involution in pair under the Poisson bracket 

mi- 

In addition, similarly to [23], we know that 

fk = ipik4'ik + '4^2k4>2k + ipskhk, l<k<N, (34) 

are integrals of motion for (j22p and (j23p . It is not difficult to verify that 3N functions {fk}k=i' 
{Fm}'^^i are in involution in pair. To show the functional independence of 3N functions {fk}k=i^ 

2N 



{Pm}m=i^ ^^6' as ™ reference [20], the technique developed by Ma et al. in 

Therefore, 3N functions {fk}k=i' {-^mlm^i are functionally independent over some region of the 
supersymmetry manifold R^^^'^^ . Now, all of the above analysis gives the following theorem. 



Theorem 1 Both the spatial and temporal constrained flows \22\l and ^23]) are Liouville inte- 
grable super Hamiltonian systems defined on the supersymmetry manifold i?^^'^^, which possess 
3N functionally independent and involutive integrals of motion {fk}k=i ^'^'^ {-Pmj^i defined by 
134^ and ( 130)) . Moreover, the formula Ii21\) provides a Bdcklund transform from the constrained 
flows (2^ and / fl3|) to the Dirac systems f77]] . 

Remark 2 The super-system on supermanifolds ^2A/|2Af Liouville integrable 126] if it pos- 
sesses M even valued conserved quantities and N odd valued conserved quantities that are inde- 
pendent and are also in involution. Furthermore, similar to the classical case, there exist a super 
analogue of Liouville' s theorem J26[ 27 J. Note that, set 9 be a odd variable in superspace, then 



fk = Gfk,k = 1,2, • • • ,N, are N odd valued conserved quantities for super-finite dimensional in 
theorem [7] because this system is involved only with GVGTl flow . 



5 Conclusions and remarks 

In this paper, we have applied the binary nonlinearization method to the super Dirac systems 
by the Bargmann symmetry constraint ()2ip . We have also shown in Theorem [1] that the nonlin- 
earized systems (I22p and ()23p are two super finite-dimensional integrable Hamiltonian systems, 
whose super Hamiltonian forms and integrals of motion have been presented explicitly. We 
would also like to emphasize that the new formula ()14p is a general result for calculating the 
variational derivative of the spectral parameter A with respect to the potential u. The crucial 
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difference between the nonlinearziation processes of the super AKNS system and the super Dirac 
system is due to the variational derivatives of A defined by formula (21) in reference [20] and 
formula (jl6p in this paper. 

We remark that the super Dirac systems ([8|) or (jlip . e.g., (jl2p . only possess super (odd and 
even) independent variables. The fully supersymmetric Dirac systems possessing both super 
dependent variables and super independent variables seem to be a very interesting object for our 
future research. For more detailed discussions on the supersymmetry theory and supersymmetric 
analysis, we would like to refer readers to reference |25j . 
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